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Optimal control modification has been developed to improve robustness to model-reference adaptive con- 
trol. For systems with linear matched uncertainty, optimal control modification adaptive law can be shown by 
a singular perturbation argument to possess an outer solution that exhibits a linear asymptotic property. Ana- 
lytical expressions of phase and time delay margins for the outer solution can be obtained. Using the gradient 
projection operator, a free design parameter of the adaptive law can be selected to satisfy stability margins. 


I. Introduction 

Various modifications were developed to increase robustness of MRAC by adding damping to the adaptive law. 
Two well-known modifications in adaptive control are the o -modification 1 and e- modification. 2 These modifications 
have been used extensively in adaptive control. Recently, a new adaptive law has been introduced that is based on 
an optimal control formulation to minimize the norm of the tracking error. 3 The optimality condition results in a 
damping term proportional to the persistent excitation. The analysis shows that the optimal control modification can 
allow fast adaptation with a large adaptive gain without causing high-frequency oscillations and can provide improved 
stability robustness while preserving the tracking performance. This paper shows that this adaptive law can exhibit an 
asymptotic linear property for systems with linear matched uncertainty using a singular perturbation argument. Due 
to the asymptotic linear property, the optimal control modification exhibits a linear mapping between the input and 
output. As a result, it is possible to compute analytically stability margins for this adaptive control scheme. With the 
use of the gradient projection operator, a free design parameter of the optimal control modification can be designed to 
satisfy stability margins for its asymptotic solution. Simulation results demonstrate the effectiveness of the method. 

II. Optimal Control Modification Adaptive Law 

Consider a linear plant with uncertainty 


x — (A T AA) x T ( B T A B) u 


(1) 


where x(t) : [0,°°) — > R" is a state vector, u{t) : [0,°°) — ■> is a control vector, A £ R" x " and B £ R nxp , p > n are 
constant and known matrices such that the pair ( A,B ) is controllable, AA and A B are constant but unknown matrices 
with an assumption that (A + AA,B + AB) is controllable, and furthermore B + AB = BA, where A £ R pxp is an 
unknown control effectiveness diagonal matrix with positive elements. 

Assumption 1: AA is a matched uncertainty whereby there exists a contant matrix Q £ R' txp such that AA = 
(B + AB) Q. . Moreover, if p > n and if rank ( (B I AB) (B + AB ) ' ] = n, then A A exists and 


qJ = (b+ab) t (b+ab) (b+ab) 1 


A A 


( 2 ) 


A nominal fixed gain controller has been designed to stabilize the plant and enable it to track a command r (t) 


Mnom 


—K x x — K r r 


Research Scientist, Intelligent Systems Division, Mail Stop 269-1, AIAA Associate Fellow 


( 3 ) 
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where r(t) : [0,°°) — + K r £ is a piecewise-continuous and bounded command vector, A — BK X £ K' iX ” is Hurwitz, 

and BK r eR" xr with r < n. 

The closed-loop nominal plant without uncertainty is 

x = A m x+B m r (4) 


This closed-loop nominal plant is then used to specify a reference model 

Xm = A m X m B m V ( 5 ) 

where x m ( t ) : [0, °°) — + R' ! is a reference state vector, and A m = A — BK X and B m = BK r . 

Since r (f ) is bounded, then x m (f ) can be shown to be uniformly bounded such that 

IM0)|| <e=> ||x m (f)|| < 5(e) , Vf > 0 (6) 

The objective is to design a full-state feedback adaptive augmentation controller to enable x(t) to follow x m (t) in 
the presence of uncertainty due to AA and A B with the following controller 


U — 11/iotn I Wad 


(7) 


where u [ld ( t ) : [0,°°) — + K p is an adaptive augmentation controller. 

Assumption 2: There exist ideal constant matrices 0* and 0* such that the following conditions are satisfied 

AA - A BK X = (B + AB) 0* T (8) 

~ABK r = (B + AB)&* J (9) 

Let ©* T = [ ©* T ©; T 1 : K px " x R pxr RP X ( m ="+ r ) and T>(x,r) = [ X T ( t ) r T (t) 1 T : K" x M r -> K m . 

Define a predictor model of the uncertain plant as 

x = A m x T B m v {B -\~ AB) 0 O -f- i^B A AB^j Uad (10) 


where x[t) : [0,°°) — > K" is an estimate of x(t) and A B ( t ) : [0,°°) — > R' !Xp is an estimate of A B. 

Choose an adaptive controller of the form 

u ad = ~@ T <t> ( 11 ) 

Then the closed-loop predictor model is described by 

i = A m x + B m r — (B + AB)0 t T>-AB© t T> + AB© t T> (12) 

where A B ( t ) = A B (t) — AS and ©(?) = © (t) — 0* are estimation errors. 

Define the tracking error dynamics as e (t) = x m ( t ) — x(t), then the tracking error equation becomes 

e = A m e + S0 t T> + AS0 T <f> + e (13) 


where B(t) = B + AB ( t ) and e (x, r) :: K” x — > K" is a second-order predictor error 

e = — AS0 t <J> 


Proposition 1: The following adaptive laws 

© = -r 0 T> [e J P - v<t> T ®B T PA m l ^j B 

AZ) t = -r B 0 T T> (e J P- v<$> T ®B J PA m 1 ^) 
is a solution that minimizes the following an infinite-time horizon cost function 

1 f tf T 

J= lim - (e — A) Q(e — A)dt 
tf^oo 2 J 0 


(14) 


(15) 

(16) 


( 17 ) 
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where A (t) : [0,°°) — > R” is a lower bound of the tracking error, r© = rl > 0 € R mxm and rg = rj > 0 £ R pxp are 
adaptive gain matrices, v > 0 € R is a free design parameter, and P = I 1 > 0 € R" x " solves 

PA m +Aj n P=~Q (18) 


where Q = Q r > 0 G R” x ". 

Proof: The cost function J is convex and represents the distance measured from a point on the trajectory of e ( t ) to 
the normal surface of a hypersphere £1= {e ( t ) G R" : ||e|| < || A|| } C S> C R". The cost function is designed to provide 
robustness by not seeking an asymptotic tracking error that tends to zero but rather a tracking error that tends to some 
lower bound away from the origin. By not requiring e (t ) — > 0 as t — > °o, the adaptation can be made to be robust. 
Therefore, the tracking performance can be traded with robustness by a suitable selection of the tuning parameter v. 

An optimal control problem can be formulated by the Pontryagin’s Minimum Principle. Define a Hamiltonian 


1 


H =2 {e ~ 


{e - A) ' 7 Q (e - A) + p T (A m e + B0 t T> + AB© T 4> + e 


where p ( t ) : [0,<=°) — > R" is an adjoint variable, then the necessary condition can be established by 

p=-VHj = -Q(e- A)- A J m p 

with the transversality condition p(tf —>■ °°) =0 since e (0) is known. 

Then the adaptive laws can be formulated by a gradient method as 3 4 

© = — r©V//gr = — r©<t>p T B 

a b t = r/jV// A/5 = f/j© 1 <pp 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 


An “approximate” solution of p (t) can be obtained using a “sweeping” method 5 by letting p = Pe + SB© <f>. Then 

d (B© t T>) 


Pe + P 


'(© T -©* T ) 


A m e + B ( © 1 -©*' ) T> + AB© T 4> + e 


+ SB© T <i> + S— - = -Q{e-A)-Aj n {pe + SB© T ^j (23) 


which yields the following equations 


P + PA w +A,' n P + Q = 0 


5 + P+A,^5 = 0 

subject to the transversality conditions P(tf ° o) =0 and S (tf — > °°) =0, ans 

_d (B© t 4>) 


QA = S- 


■ —P{B + AB) ©* T <J> 


dt 


Introducing a time-to-go variable x = tf — t, then the Lyapunov differential equation (24) becomes 

dP T 

— — — t -PA m +A n P + Q = 0 
d x 


(24) 

(25) 


(26) 


(27) 


subject to the initial conditions in time-to-go P( 0) = 0, for which the existence and uniqueness of its solution is 
well-established. 

. It follows that Eq. (25) also has a stable, unique solution in time-to-go since A m is Hurwitz 


-f T+ PB + AlS = 0 

subject to S (0) = 0. 

The solutions of P (t) and S ( x) as tf — * °° approach their steady-state constant solutions 

PA m +Aj n P = — Q 
S = -atJp 


(28) 


(29) 

(30) 
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Without any loss of generality, a free design parameter v > 0 € K. is introduced as a gain to allow for adjustments 
of the modification term in the adaptive law, where v = 1 corresponds to an optimal solution. Thus 

S = —vA m J P (31) 

The adjoint p is then obtained as 

p = Pe- vA m J PB® J <t> (32) 

Substituting p into the gradient adaptive laws yields 

© = -r 0 T> (e T P- v<t> T ®B T PA m 1 ^) B (33) 

AZF = -r B ® T T> (e T P - v<t> T ®B T PA m l ^J (34) 

Theorem 1: The adaptive laws (15) and (16) result in stable and uniformly ultimately bounded tracking error e ( t ) 
for all (e (0) , © (0) , A B (0)) € B a with an ultimate bound 


P = 


where 


I hnax (P) r 2 + hnax (Tq*) K 2 + Xmax (T^ 1 ) V~ 
^ min (P) 

r= 2 H^li £q 

hum ((2) 

2||PA m 1 ||0 o 
^ min (^m t qa 

m 

2\\PA- l \\AB 0 
^ min i^m t QA m 


K = 


V = 


with £ 0 = supjg^gur || e (x,r)||, ©o = ||©*||, and A B 0 = ||AB||. 
Proof: Choose a Lyapunov candidate function 


V = e T Pe + trace 


^©^g, 1 ©^ + trace 


ABF b 1 AB t 


Evaluating V yields 


(35) 


(36) 

(37) 

(38) 


(39) 


V = e T (A,„P + PA m ) e + 2e T P ^B© T 4> + AB© T <t> + £^ -2trace © T 4> ( 

— 2trace 

Using the trace identity trace ( A r B ) = BA T , V can be written as 


© T <f> ( e T P- v<t> J ®B J PA7J ) B 


A B® T <t> ( e J P — v<l> 1 ®B 1 PA 


T c\£T da — 1 


(40) 


V = ~e J Qe + 2e J P (z?© T 4> + AB© T <1> + fi) -2 e T PB© T <f> 

+ 2v<t> T ©B T PA," 1 fi© T <t> - 2e T PAB® T <P + 2v<S> T ®B T PA m 1 AB© T 4> (41) 


The sign-definiteness of the term PA m 1 is now considered. Recall that a general real matrix G is positive (negative) 
definite if and only if its symmetric part M = ^ (G + G T ) is also positive (negative) definite. Then, by pre- and post- 
multiplication of Eq. (18) by A ln and ,4 " 1 , respectively, PA" 1 can be decomposed into a symmetric part M and 
anti-symmetric part N as 


PA m l =M + N 


(42) 


where 

M =\ (A-TP + PA" 1 ) = - ] -A m J QA-' 
N= l -(PA- 1 -A~ T P). 


(43) 

(44) 
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Since the symmetric part M < 0, then PA m 1 < 0. Thus, V becomes 
V = -e J Qe+2e J Pe + 2v<P J (©* + ©) B 1 (M + N) B® r <i> + 2v<P T ® (AB + AB) J (M + N) AB® t <P (45) 

Using the property y T Ny = 0 for an anti-symmetric matrix N, V is reduced to 

V = -e T Qe + 2e T Pe + 2v<t> T ®*B T PA m 1 B@ T <i> - v<t> T ®B J A m r QA m l B® T <t> 

+ 2v4> T 0AB T PA m ' AB© t 4> - v^ t @AB t A- j QA- 1 AB® t ^ (46) 

which is bounded by 


v < - Ikll [Kin (Q) Ikll -2 \\P\\ eo] - v ||£|| ||0|| 2 ||© 

- V 


Kin Q^m 


© -2 PA' 1 ©o 


2 IIo|| 2 ||ab| 


Kin {^m Q^m 


AB -2 PA- 1 ABo 


Let 

where 


B r = {(e,©,AZ?) el*x R mxp x R nxp : ||e|| < r or 1 1 © 1 1 < ifor ||AB|| < u} 

2 \\P\\ £q 


Ikll < r = 

l©ll <K = 


Afi < v = 


^ min ( Q) 

2||M m 1 ||0 o 
tin QAm 
2\\PA m l \\AB Q 
Klin (\ B^ A m QA m 


B 


(47) 

(48) 

(49) 

(50) 

(51) 


It follows that V < 0 for all (e,©,A B) £ Br — B r , where Br = {e £ K” : ||e|| < R} C Let Bp be the smallest 
subset that encloses B r , then there exists j3 > 0 where 


P — Kax ( P)r 2 + hnax (T©*) K 2 + Kax ') V~ 


such that 


B r cBp = {(<?,©, A B) £R n xR mxp xR nxp :V < p} 

Let B a be the largest subset enclosed by Br, then since ||e|| < R in Br, there exists a > 0 where 

Kin ( P ) Ikll 2 < Kin ( P ) Ikll 2 +Kin (r© 1 ) ||©|| 2 + Kin (r B 1 ) ||AB|| 2 < V < Kin ( P)R 2 = a 

such that 

B a = { (e,©, AB) £R n xR mxp xR" xp :V < a} CB R 

Then for a solution to be uniformly bounded, the set containment is as follows: 


(52) 

(53) 

(54) 

(55) 


B r CBp C B a C Br 


This implies 
Therefore 


P < a Kax ( P ) r 2 + Kax (r © 1 ) K 2 + Kax (r^ U 2 < Kin (P) R 2 


R> 


I ^max (-P) T 2 + Xmax (^@ 1 ) K* 2 + X max (T B 1 ) U 2 
Kim ( P ) 


= P 


where p is the smallest value of R. 

Then p is the ultimate bound of e ( t ) such that 


(56) 

(57) 

(58) 


r<|kll<P<^ (59) 

Since V < 0 for all ( e , © . AB ) £ Br — B r , therefore V is a decreasing function of time outside of B r . Thus, if 
(e (0) ,0(0) ,AB(0)) £ B a , then according to Theorem 5.1 of Ref., 6 the solution will eventually enters Bp after a 
finite time t = T (independent of (e (0) ,©(0) , AS(0)) and a) and remain inside for all t > T. Therefore, e (t) is 
uniformly ultimately bounded with an ultimate bound p . 
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III. Asymptotic Linearity 


Note that the tracking error equation can also be written as 

e=A in e-(B + AB)@* J <t> + B® T <t> (60) 

Now consider a time-scale separation problem when ||F©|| = ® ||A m || and ||F ^|| = y ||A„, || where e is 

a small parameter. Then © (t) and A B (t) are said to be a fast state and e ( t ) is a slow state. To decouple the fast 
and slow states, a time-scale separation can be performed by applying the singular perturbation method using a time 
transformation 

X = £t (61) 

where T is a slow time variable. 

Then, the system is transformed into a singularly perturbed system as 

— = - \A m e - (B + A B) ©* T $ + B© T $| (62) 

d T e L J 

e— ^-—<p(e T P-v<i> T ®B T PA m l )B (63) 

dx e V ) 

e ^K-^® T *(e T P-v* T ®# T /A-*) (64) 

The Tikhonov’s theorem can be used to approximate the solution of the singularly perturbed system with the 
solution of a “reduced-order” system by letting e — > 0. 7 Then, 0 (e. e) and A B are on a fast manifold. Thus, the 
reduced-order system is given by 

£ 2 ^ « -70$ (e J P- v<P T ®B T PA m 1 ^) B => (^© T $) l0) = ^p- l A m l Pe {0) (65) 

e (0) = A,„e (0) - (B + A B) ©* T $ (0) + (b@ t $) ( ° } (66) 

where (t) and (B© $) ^ (t) are the “outer” solutions of the singularly perturbed system. 

The term “outer” is in connection with the concept of “inner” or “boundary layer” and “outer” solutions which 
have the origin in boundary layer theory due to Prandtl. The “inner” or “boundary layer” solution for this system is 
obtained from 

@(0 = -r©$® (e (i)T P-V ($ t ©B t ) PA m 1 j B ( ') (67) 

A& 1 = -r b ©W t $W (eM J P- v (<P T ®B j y ] PA~ l ^J (68) 

e^L = A m e — (B + AB) ©* T <t> + B©$ => A m e^ = (B + AB) ®* T - (b@ t (69) 

The general solution of the system is then expressed as 

e (; t ) = e (0) (t) + e (l) (t) - e M AE (t) (70) 

where cmae (?) is a correction term by a matched asymptotic expansion method applied to both the inner and outer 
solutions. 8 The outer solution is in fact the asymptotic solution of the original system as t — > For the singularly 
perturbed system comprising e !l> > (t) and (B© <J>) '° (t), the outer solution of the tracking error is then determined by 

e (0) = A m + ^P~ l Aj t P+(AA — ABK x ) e (0) - (AA - ABK X ) x m + ABK r r (71) 

which interestingly enough is an LTI equation. 

Since the asymptotic tracking error equation for the outer solution is LTI, many standard LTI analysis tools can 
be used to estimate the asymptotic behavior of the original system. In particular, the free design parameter v can be 
selected to provide a notion of stability margin for the outer solution. 
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The projection operator can be applied to the adaptive laws with a priori knowledge of the bounds on A A and A B. 
Using these a priori bounds, the free design parameter v can be selected to guarantee stability robustness. Then the 
adaptive laws are modified by the projection operator as follows: 

© = -r 0 Proj (©,4> (e T P- v4> T ©B T PA m 1 ) B) (72) 

AS t = -r B Proj (AB,© t <I> (e J P- v4> T ©B T PA m 1 ^ (73) 

The free design parameter v can now be estimated to give the asymptotic closed-loop tracking error a desired 
measure of stability robustness such as phase or time delay margin. One method for estimating a MIMO phase or time 
delay margin is based on a matrix measure method 9 which is as follows 
In the presence of an input delay, the closed-loop tracking error becomes 


e (0) (t) = (A + AA)e {0) (t) - 


(B + AB)K X - ^P~ l Aj,P 


e (0) (: t - t d ) - (AA — A BK X ) x m (t-t d )+ A BK r r (f - t d ) (74) 


The MIMO phase and time margins can be estimated by 


A 1 A (A + AA) + A (; [(B + AB) K x - l -P~ x Aj n P ] ) 

(2) — COS « 

\\{B + AB)K X -\P-'AIP\ 


~p(-jA- jAA)+\\(B + AB) K x - ^P- l Aj,P\ \ 
where fl as the maximum eigenvalue of a symmetric part of a general complex matrix C 


A (C) = K 


ax 


i^f) 


(75) 

(76) 


(77) 


where C* is the complex conjugate of C. 

Utilizing this equation, the free design parameter v can then be estimated to provide a measure of a phase margin 
in an asymptotic sense. 

Example: Instability of adaptive systems with unmodeled dynamics is a well-known phenomenon. 10 Consider 
a stable first-order plant coupled with a second-order unmodeled plant with two highly damped poles which could 
represent a structural mode 

x = — x + 2u — O.ly 
y + 2£(» n y + (Ony = 7x 

where £ = 5 and co„ =10. 

The reference model is 

x m = — 2x m + 2 r 

where the reference command is a step input with a small sinusoidal variation at the same frequency as G)„ and a small 
noise signal 

r = 1 + 0. 1 sin c o n t 


The controller is given by 


where <t> (t) 


x{t) 


r{t) 


-i T 

and © (t) 


u = © T <t> 


0,(0 


e r (t) 


-I T 

is updated by 


©(f) =r(<I><?-v<M> T ©) 

with e(t) =x ( t ) — x m (t), © (0) = 0, and T = I. 

The conventional MRAC, i.e., v = 0, leads to an instability as shown in Fig. 1 even though the open-loop plant is 
stable. A value of v = 80 is chosen, the optimal control modification is able to produce a stabilizing controller, as seen 
in Fig. 2. Thus, this illustrates the robustness of the optimal control modification. As a comparison, the e-modification 
is used instead for the update law 

© = r(<t>e — jU |e|0) 
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with = 80. The adaptive signal is stable as shown in Fig. 3. 
The asymptotic tracking error equation for fast adaptation is 


e 

u — > — 

v 


The asymptotic closed-loop plant can be expressed as 


2 0.7 \ 4 r 

v + s 2 + 2^co n s + (o}j ' ' v(s + 2) 


As t — > oo 

x 2 a )}, 

r (v-2)fi)2 + 0.7v 



Fig. 1 - Instability due to MRAC, V = 0 



Fig. 2 - Stable Adaptation with Optimal Control Modification, v = 
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Fig. 3 - Stable Adaptation with e-Modification, fj. = 80 

Figure 4 illustrates the linear mapping between x and r for three different multipliers 0.1, 0.5, and 1.0 that scale 
the input signal. The predicted asymptotic ratio of x to r is -0.0255, which agrees surprising well with the simulation 
results. In contrast, Fig. 5 shows that the e-modification results in a nonlinear mapping between x and r. 



Fig. 4 - Input-Output Linear Mapping by Optimal Control Modification 
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Fig. 5 - Input-Output Nonlinear Mapping by e- Modification 


IV. Flight Control Simulation 

Consider a longitudinal pitch dynamical model of an aircraft 


mV 


CmqSc 


2V 


0 

n a ( ‘ 

2V 


Cm^qSS 2 


0 

0 


a 


1 

0 


6 

= 

0 

lyy . 


4 _ 



CrnqSc 


mgj — C La qS —mgy mV-—^ 

0 0 1 


C, 


m a 


0 


C mq qSc 2 
2V 



a 


i 

<o 


e 

+ 

0 


q . 


1 

i 


(78) 


A numerical model for a full-scale generic transport model (GTM) at Mach 0.8 and 30,000 ft with the flight path 
angle y = 0 is given by 


a 


’ -0.7018 

0 

0.9761 


a 


' -0.0573 " 

9 

= 

0 

0 

1 


e 

+ 

0 

q . 


-2.6923 

0 

-0.7322 


. q . 


-3.5352 


A desired reference model of the pitch attitude is given by 

0m + 2 C©,,0, 


(o},d m = mlr 


(79) 


where £ = 0.85 and ft),, = 1.5 rad/sec are chosen to give a desired handling characteristic. 

T 


Let x = 


a 0 q 


*3 


«31 2£o„ + (333 


and it = 8,.. A nominal controller is designed as u nom = —K x x — k r r where K x = 

and k r = — (0„ = 0.6365. The closed-loop eigen- 


0.7616 -0.6365 -0.5142 
values are —0.6582 and —1.2750 ±0.7902/. The nominal closed-loop plant is then chosen to be the reference model 
as 


a 


' -0.6582 

-0.0365 

0.9466 


a 


' 0.0162 " 

9 

= 

0 

0 

1 


e 

+ 

0 

q . 

m 

0 

-2.2500 

-2.5500 


q . 

m 

1 
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Suppose an in-state uncertainty and in-control uncertainty exist in the form of 



0 


-0.0573 





0 

0 

0.0119 

A B = 

0 

, AA = 

0 


0 0 

-0.2071 

= 

0 

0 

0 


1.7676 


-1.7676 





0 

0 

0.3661 



s. 










B+AB 


which represents a 50% reduced pitch damping and control effectiveness and 11% increase in ^-contribution to a. 

Figures 6 is a plot of estimates of phase and time delay margins from the asymptotic solution computed from Eqs. 
(75) and (76) as a function of v. Note that the phase margin with AA and AS decreases steadily as v increases, while 
the time delay margin reaches a maximum at about v = 1. Thus, for practical design purposes, v should be kept 
between 0 and 1 . A large value of v produces a better time delay margin, but also results in a poorer steady-state 
tracking. 

Suppose there exists a 0.02 sec time delay at the input due to a latency in the flight control system. Then from Fig. 
6, a value of v = 0.1 is selected to provide a time delay margin of about 0.027 sec. It should be noted that the time 
delay margin as computed from Eq. (76) is a conservative estimate corresponding to a large adaptive gain value. So 
the actual time delay margin should be greater than the estimated value. 



Fig. 6 - Phase and Time Delay Margin Estimates 

Figures 7, 8, and 9 show the a, 0, and q responses due to the baseline controller. With no adaptation, the plant 
is not able to track the reference model very well. For the standard MRAC, an adaptive gain of T = 1000/ is used to 
enable a tight tracking of the pitch attitude reference signal. The control signal is highly oscillatory during the first 15 
sec as evidenced in Fig. 10. The simulation result of the time delay margin of the system with the standard MRAC 
is observed to be 0.07 sec. For the optimal control modification, the design parameter v = 0.1 is used with the same 
adaptive gain. The response was more improved without high frequency signals. The observed time delay margin 
of the optimal control modification adaptive law is 0.3 sec which is significantly greater than that with the standard 
MRAC. 
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Fig. 7 - a, 0, and q Responses with Baseline Controller 




Fig. 8 - a, 9, and q Responses with Standard MRAC 
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Fig. 9 - a, 0, and q Responses with Optimal Control Modification (OCM) 



Fig. 10 - Elevator Deflection 


V. Conclusions 

This study shows that the recently developed optimal control modification exhibits a linear asymptotic property 
for systems with linear matched uncertainty. Using a singular perturbation argument, the outer solution of the adaptive 
system is governed by an asymptotic linear time invariant system. As a result, the notions of phase and time delay 
margins for this asymptotic system can be established analytically. With the use of the projection operator, the optimal 
control modification can be designed to satisfy stability margins in a linear sense by a suitable selection of the free 
design parameter of the adaptive law. This approach might provide a path toward easing technical challenges with 
verification and validation of adaptive control systems. 
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